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Similarly to quantum states, also quantum measurements can be "mixed", corresponding to a 
random choice within an ensemble of measuring apparatuses. Such mixing is equivalent to a sort of 
hidden variable, which produces a noise of purely classical nature. It is then natural to ask which 
apparatuses are indecomposable, i. e. do not correspond to any random choice of apparatuses. This 
problem is interesting not only for foundations, but also for applications, since most optimization 
strategies give optimal apparatuses that are indecomposable. 

Mathematically the problem is posed describing each measuring apparatus by a positive operator- 
valued measure (POVM) , which gives the statistics of the outcomes for any input state. The POVM's 
form a convex set, and in this language the indecomposable apparatuses are represented by extremal 
points — the analogous of "pure states" in the convex set of states. Differently from the case of states, 
however, indecomposable POVM's are not necessarily rank-one, e. g. von Neumann measurements. 

In this paper we give a complete classification of indecomposable apparatuses (for discrete spec- 
trum), by providing different necessary and sufficient conditions for extremality of POVM's, along 
with a simple general algorithm for the decomposition of a POVM into extremals. As an interesting 
application, "informationally complete" measurements are analyzed in this respect. The convex set 
of POVM's is fully characterized by determining its border in terms of simple algebraic properties 
of the corresponding POVM's. 
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I. INTRODUCTION 



Measurements are the essence of any experimental sci- 
ence. At extreme sensitivities and precisions they become 
the true core of quantum mechanics. For any practical 
need, a measurement can be always regarded as the re- 
trieval of information about the "state" of the measured 
system. However, due to the no-cloning theorem]^, ||, 
even for an elementary system — e.g. a single harmonic 
oscillator or a spin — it is impossible to recover a com- 
plete knowledge of the state of the system from a single 
measurement H without prior knowledge. Then, since 
in quantum mechanics different incompatible measure- 
ments can be performed in principle, one is faced with 
the problem of which measurement should be adopted 
for accomplishing a specific task, and which strategy of 
repeated measurements would be the most statistically 
efficient. These are the basic issues of the operational 
viewpoint of Quantum Estimation Theory Q. 

A measurement on a quantum system ||^ returns a 
random result e from a set of possible outcomes E = {e = 
1, . . . N}, with probability distribution p{e\p) depending 
on the state p of the system in a way which is distinctive 
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of the measuring apparatus according to the Born rule 

p(e|p) = Tr[pPe]. (1) 

In Eq. (^ Pe denote positive operators on the Hilbert 
space H of the system, representing our knowledge of 
the measuring apparatus from which we infer informa- 
tion on the state p from the probability distribution 
p{e\p). Positivity of Pg is needed for positivity of p(e|p), 
whereas normalization is guaranteed by the completeness 
^ggg Pe = /. In the present paper we will only consider 
the simple case of finite discrete set E. More generally, 
on has an infinite probability space E (generally contin- 
uous), and in this context the set of positive operators 
{Pe} becomes actually a positive operator valued mea- 
sure (POVM). Every apparatus is described by a POVM, 
and, reversely, every POVM can be realized in principle 
by an apparatus |, |, |, |, |, 

The linearity of the Born rule (|1|) in both arguments p 
and Pe is consistent with the intrinsically statistical na- 
ture of the measurement, in which our partial knowledge 
of both the system and the apparatus reflects in "con- 
vex" structures for both states and POVM's. This means 
that not only states, but also POVM's can be "mixed", 
namely there are POVM's that give probability distri- 
butions p{e\p) that are equivalent to randomly choosing 
among different apparatuses. Notice that mixed POVM's 
can also correspond to a single measuring apparatus, not 
only when the apparatus itself is prepared in a mixed 
state, but also for pure preparation, as a result of discard- 
ing (tracing out) the apparatus after a unitary interaction 
with the system. Clearly, such mixing is itself a source of 
"classical" noise, which can be in principle removed by 
adopting an indecomposable apparatus in the ensemble 
corresponding to the mixed POVM. It is then natural to 
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ask which apparatuses are indecomposable, i. e. "pure" 
in the above sense, or, mathematically, which POVM's 
correspond to extremal points of the convex set. The 
classification of such apparatuses is certainly very use- 
ful in applications, since most optimization strategies in 
quantum estimation theory j|, |^ correspond to minimize 
a concave (actually linear) function on the POVM's con- 
vex set — so-called "cost- function" — whence leading to an 
optimal POVM which is extremal. 

Surprisingly, extremality for POVM's generally doesn't 
mean to be rank-one, as for "pure" states. In other 
words, indecomposable POVM's are not necessarily re- 
alized by von Neumann measurements. Indeed, as we 
will see in this paper, there are rank-one POVM's that 
are not extremal, whereas, on the opposite side, there 
are higher-rank POVM's which are extremal. Moreover, 
whenever the optimization problems have additional lin- 
ear constraints — e. g. for covariant POVM's, or for 
fixed probability distribution on a given state — the cor- 
responding subset of POVM's is a lower dimensional con- 
vex set corresponding to a section by a hyperplane of the 
complete POVM's convex set, with boundary equal to 
the section of the original boundary, and whence with 
extremal points that belong to the boundary of the con- 
vex set of all POVM's. For this reason, also the boundary 
of the POVM's convex set is interesting in practice, since 
POVM's that are optimal (for a concave cost function) 
with an additional linear constraint generally are non ex- 
tremal, but still belong to the boundary. One can also 
argue that POVM's which lie inside faces of the convex, 
physically exhibit a different degree of "classical" noise 
in relation with the dimensionality of the face. Notice 
that one should not imagine the POVM's set as a poly- 
tope, since, on the contrary the set is "strongly convex" , 
namely the extremal points are not isolated, but lie on 
continuous manifolds. 

In the present paper we address the problem of appa- 
ratus dccomposability along three lines of attack: a) by 
providing simple necessary and sufhcient conditions for 
extremality of POVM's; b) by establishing the complete 
structure of the POVM's convex set via the character- 
ization of its border in terms of algebraic properties of 
the POVM's; c) by providing a simple general algorithm 
for the decomposition of POVM's into extremals. For 
simplicity, the whole paper is restricted to the case of 
discrete spectrum. In Section |l|, after clarifying the gen- 
eral features of convex combinations of POVM's, using 
the method of perturbations we derive three different if- 
and-only-if conditions for extremality, along with some 
corollaries giving easy useful conditions, only necessary 
or sufficient, that will be used in the following. Section 
III exemplifies the results of Section || in the case of a sin- 
gle qubit. Section!^ presents the characterization of the 
border of the convex set in terms of algebraic properties 
of POVM's. Section shows that for every dimension 
d — dim(H) there is always an extremal POVM with 
maximal number TV of elements N = d^ , correspond- 
ing to a so-called informationally complete P0VM[|l2j. 



After summarizing the results in the concluding Section 
VI , Appendix ^ reports the algorithm for decomposing 
a POVM into extremals. 



II. CONVEXITY AND EXTREMALITY OF 
POVM'S 

Let's denote by the convex set of POVM's on a 
finite dimensional Hilbert space H, with a number N of 
outcomes E = {l,...,Af}. We will represent a POVM 
in the set as the vector P G P = {Pi,...,P/v} 

of the N positive operators Pf.. The fact that the set 
■!^Ar is convex means that it is closed under convex lin- 
ear combinations, namely for any P',P" G also 
P = pP' + (1 - p)P" e with < p < 1 — i.e. p is 
a probability. Then, P can be also equivalently achieved 
by randomly choosing between two different apparatuses 
corresponding to P' and P", respectively, with probabil- 
ity p and 1 — p, since, the overall statistical distribution 
p{e\p) will be the convex combination of the statistics 
coming from P' and P". Notice that contains also 
the set of POVM's with a strictly smaller number of out- 
comes, i. e. with E' C E. For such POVM's the elements 
corresponding to outcomes in E \ E' will be zero, corre- 
sponding to zero probability of occurrence for all states. 
Clearly, for TV < M one has Q C ^, where by 
3^ we denote the convex set of all POVM's with any (gen- 
erally infinite) discrete spectrum. The extremal points 
oi represent the "indecomposable" measurements, 
which cannot be achieved by mixing different measure- 
ments. Obviously, a POVM which is extremal for is 
also extremal for 3^m with M > N, whence it is actu- 
ally extremal for and we will simply name it extremal 
without further specifications. 

Let's start with a simple example. Consider the fol- 
lowing two-outcome POVM for a qubit 

P=(i|0)(0|, i|0)(0| + |l)(l|). (2) 

By defining D = i(|0)(0|, -|0)(0|), the two vectors P± = 
P ± D correspond to the following different POVM's 

P+ = (|0)(0|, |1)(1|), P--(0,/), (3) 

which, intuitively are extremal, whereas P is not, since 
P = ^P+ -|- ^P- (we will see in the following that P± 
are indeed extremal). 

Now the problem is how to assess when a POVM is 
extremal. Looking at the above example, one notices that 
the non extremality of P is equivalent to the existence 
of a vector of operators D ^ such that P± = P ± D 
are POVM's, because in this case P can be written as a 
convex combination of P-|_ and P_. The non existence 
of such vector of operators D 7^ is also a necessary 
condition for extremality of P, since for a non extremal 
P there exist two POVM's Pi ^ P2 such that P = 
iPi -I- iP2, whence D = Pi - P2 7^ 0. 
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This leads to the method of perturbations for estab- 
Ushing extremahty of a point in a convex set, which in 
the present context will be the following. 

Method of perturbations. We call a nonvanishing D a 
perturbation for the POVM P if there exists an e > 
such that P ± eD arc both POVM's. Then a POVM is 
extremal if and only if it doesn't admit perturbations. 

From the definition it follows that perturbations for 
POVM's are represented by vectors D of Hermitian op- 
erators (for positivity of P ± eD) and with zero sum 
X^eeE^e = (for normalization of P ± eD). Specifi- 
cally, D is a perturbation for P if for some e > one 
has Pe ± eDe > Ve G E. Notice that the condition 
Pe±eDe > is equivalent to e\De\ < Pe, and a necessary 
condition for e\De\ < Pe is that Ker(Pe) C Ker(_De), or 
equivalently Supp(De) C Supp(Pe) [the support Supp(X) 
of an operator X is defined as the orthogonal comple- 
ment of its kernel Ker(X)]. In fact Ker(_De) = Ker(|Z?e|), 
and for a vector lip) £ Ker(Pe) with \-4>) ^ Ker(|Z?e|), 
one would have (■(/'|(Pe - e|£'e|)|V') = -£("011-06 11-0) < 0, 
contradicting the hypothesis. Clearly, the Hermitian op- 
erators De can be taken simply as linearly dependent — 
instead of having zero sum — i.e. \eDe = for non- 
vanishing Ae, and, moreover, one can consider more gen- 
erally complex operators De with Supp(-De) U Rng(_De) C 
Supp(Pe) Ve S E, and satisfying e\De\ < Pe- In fact 
Pe ± eDg > is satisfied Ve G E by the set of Hermitian 
operators D'^ = XeDe + X^Dl, for which J2e ^'e = 0- 

The above considerations show that what really mat- 
ters in assessing the extremahty of the POVM P = {Pe} 
is just a condition on the supports Supp(Pe), correspond- 
ing to the following theorem. 

Theorem 1 The extremality of the POVMP is equiva- 
lent to the nonexistence of non trivial solutions D of the 
equation 

J^De^O, Supp(7?e) U Rng(i^e) C Supp(Pe) Ve G E. 

eeE 

(4) 

The above condition can be made explicit on the Pg 
eigenvectors corresponding to nonzero eigen- 

value, which therefore span Supp(Pe). Then, Eq. (^) 
becomes the linear homogeneous system of equations in 

the variables Dnm = {vn''\De\vm) 

rank(_Pe) 

eGE nm — 1 ^ 

^ = Vn, m, Ve e E, 

namely the following version of Theorem |^ 

Theorem 2 A POVM P = {PejeeE is extremal iff the 
operators \vn''){vm^\ made with the eigenvectors of Pe are 
linearly independent Ve G E and Vn, m = 1, . . . , rank(Pe) . 



Theorem || is the characterization of extremal POVM's 
given by Parthasaraty in Ref. ||l^ in a C*-algebraic set- 
ting. Notice that instead of the eigenvectors {Ifi'^^)} one 
can more generally consider a non orthonormal basis, 
which can be useful in numerical algorithms. 

Another interesting way to state Theorem |l] is in terms 
of weak independence of orthogonal projectors Ze on 
Supp(Pe), where we call a generic set of orthogonal pro- 
jections {Ze}eeE weakly independent if for any set of 
operators {PejeeE on H one has 

E ZeTeZe = ^ ZeTeZe = 0, Ve G E. (6) 
e6E 

Notice that since the orthogonal projectors are in one-to- 
one correspondence with their supporting spaces, the no- 
tion of weak independence can be equivalently attached 
to the supporting spaces Supp(Pe). This leads us to the 
alternative characterization of extremal POVM's 

Theorem 3 A POVM P = {PejeGE is extremal iff the 
supports Supp(Pe) are weakly independent for all e G E. 

The proof is straightforward if one considers that any 
perturbation D for P can be written as De = ZeDeZe 
with the constraint 

E ZeDeZe = . (7) 

e 

Then the theorem simply says that the only allowed per- 
turbation for an extremal POVM is the trivial one. 

Some corollaries relevant for applications follow imme- 
diately from the main Theorem |l| or its equivalent ver- 
sions. 

Corollary 1 Por ^^^^ dim[Supp(Pe)]2 > the POVM 
P is not extremal. 

This means that a POVM with more than nonvanish- 
ing elements is always decomposable into POVM's with 
less than d^ elements. For the case of elements Theo- 
rem |l| also implies that 

Corollary 2 An extremal POVM with d'^ outcomes must 
be necessarily rank- one. 

In fact, for larger rank we would have more than d^ eigen- 
vectors \vn^), and the operators \vn^){vm \ cannot be lin- 
early independent. From Theorem m it also follows that 
if some elements of the POVM P have non-disjoint sup- 
ports, then P is not extremal, or, equivalently 

Corollary 3 An extremal POVM P = {Pe} must have 
all supports Supp(Pe) mutually disjoint. 

Precisely, we call two linear spaces A and B disjoint when 
A n B = {0} is the null vector. It is worth noticing 
that two linear spaces that are disjoint are not necessar- 
ily orthogonal, whereas, reversely, two orthogonal spaces 
are clearly disjoint. We emphasize that the condition of 
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Corollary ^ is only necessary. Indeed, it is easy to envis- 
age a POVM that satisfies the above condition without 
being extremal, e. g. a rank-one POVM for d = 2 with 
five elements corresponding to the vertices of a pentagon 
in the Bloch sphere (see Section III for extremal POVM's 
for qubits). 

Other obvious consequences of Theorem || are the fol- 
lowing 

Corollary 4 Orthogonal POVM's are extremal. 

Corollary 5 A rank-one POVM is extremal if and only 
if its elements Pe are linearly independent. 

Corollary ^ states that for dimension d an extremal 
POVM can have at most d^ non-null elements (in Sec- 
tion ^ we will show that such an extremal POVM always 
exists). Here we can immediately conclude that 

Corollary 6 A POVM with d^ elements is necessarily a 
rank-one "informationally complete" POVM. 

By definition, an informationally complete POVM P = 
{Pe} has elements Pe which span the space of all 
operators on H, thus allowing the estimation of any en- 
semble average using the same fixed apparatus. The fact 
that a POVM with d^ elements necessarily is rank-one is 
stated in Corollary whereas Corollary ^ assures that 
all POVM's elements are linearly independent, whence 
the set of d^ linearly independent operators Pe is obvi- 
ously complete for dimension d. In Section ^ we will give 
an explicit example of such an extremal informationally 
complete POVM. 

For rank greater than one we have only the necessary 
condition 

Corollary 7 // the POVM is extremal, then its non- 
vanishing elements are necessarily linearly independent. 

In fact, according to Theorem || the projectors on the 
eigenvectors must be linearly independent, whence also 
the operators Pe- Indeed, for linearly dependent ele- 
ments there exist coefficients Ae not all vanishing such 
that X^eGE ''^e-^e = 0, and without loss of generality we 
can take -1 < A. < 1. Therefore, the POVM can be 



written as convex combination P — +^P+, with 



P^ = (1 ± Xe)Pe, and ^ P+ (since the Ae's are not 
all vanishing). 

In appendix ^ we report an algorithm for decomposing 
a given POVM into extremal. 



III. EXTREMAL POVM'S FOR QUBITS 

Using the above results we will give a classification 
of extremal POVM's for qubits. In this case. Corollary 
|l| implies that the extremal POVM's cannot have more 
than 4 elements, and that, apart from the trivial POVM 
P = /, they must be made of rank-one projectors (other- 
wise Supp(Pe) for different e G E would not be mutually 



disjoint). Now, upon writing the POVM elements in the 
Bloch form 



Pe ^ ae{I + fie ■ S) . 



(8) 



the constraints for normalization and positivity read 



tte > , 



ae 



(9) 



The case of two outcomes corresponds simply to the usual 
observable Pe = |e),-j fi(e|, e = 0, 1, with |e)fj eigenvector 
oi fi ■ a corresponding to the eigenvalues -1-1, —1, respec- 
tively. In fact, for two outcomes one has aofio+aifii — 0, 
namely = —fii n, and necessarily ao ai = ^. We 
now consider the case of 3 and 4 elements. By Theorem 
a necessary and sufficient condition for extremality is 



Jeae{I -I- rie • (?) = 7e = Vc G E, 



(10) 



eeE 



or, equivalently, 



^leae 

eeE 



0, JeCtene 







7e = 0VeeE. (11) 



For three outcomes Eq. (||) implies that {aefieleeE rep- 
resent the edges of a triangle, and thus the second condi- 
tion in Eq. (|ll| ) is satisfied iff 7e = 7 is independent of e. 
Then the first condition is satisfied iff 7 = 0. Therefore, 
all three outcomes rank-one POVM's with pairwise non 
proportional elements are extremal. For four outcomes 
we can see that for an extremal POVM the correspond- 
ing unit vectors {ne}e<£E cannot lie on a common plane. 
Indeed, divide the four vectors {aene}eeE into two cou- 
ples, which identify two intersecting planes. Then, the 
third condition in Eq. (^ implies that the sums of the 
couples lie on the intersection of the planes and have 
the same length and opposite direction. If we multiply 
by independent scalars 7e the two elements of a couple, 
their sum changes direction and lies no more in the in- 
tersection of the two planes, and the second condition 
in Eq. (^Tj) cannot be satisfied. Therefore, the two ele- 
ments of the same couple must be multiplied by the same 
scalar, which then just rescales their sum. Now, when the 
rescaling factors are different for the two couples the two 
partial sums don't sum to the null vector anymore. Then 
necessarily 7e = 7 independently of e. In order to satisfy 
also the first condition in Eq. (^l]) we must have 7 = 0. 
On the other hand, if the four unit vectors lie in the 
same plane, a non-trivial linear combination can always 
be found that equals the null operator, hence the POVM 
is not extremal. By the first condition in Eq. ( pj] ) we 
have indeed 



70 ao 



E 

e=l,2,3 



Jeae 



then the second condition can be written 

y^ 7eae(ne - ?1o) = . 
e=l,2,3 



(12) 



(13) 
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Now, either we have a couple of equal vectors fie, or the 
three vectors rie — in a two dimensional plane are lin- 
early dependent. However, in both cases the POVM is 
not extremal, because in the former case two elements 
are proportional, while in the latter a non trivial triple 
of coefficients 7e satisfying Eq. ( |l3| ) exists. 

Notice that, the three and four-outcomes POVM's are 
necessarily unsharp, i.e. there is no state with probability 
distribution p{e\p) = de,e for a fixed e. They provide ex- 
amples of un-sharp POVM's with purely intrinsical quan- 
tum noise. 



IV. THE BOUNDARY OF THE CONVEX SET 
OF POVM'S 

In this section we generalize the results about ex- 
tremality, and give a full characterization of the elements 
on the boundary of the convex set of n-outcomes POVM's 
on the Hilbert space H. Let start from an intuitive geo- 
metrical definition of the boundary of a convex set. Con- 
sider for example a point lying on some face of a polyhe- 
dron. Then there exists a direction (e. g. normal to the 
face) such that any shift of the point along that direction 
will bring it inside the convex set, while in the opposite 
direction it will bring the point outside of the convex. 
In mathematical terms, consider a convex set and an 
element p e Then, p belongs to the boundary of 
'j^ if and only if there exists such that 

p + e(g-p)e<r, p^e{q~v)^'^. Ve G [0, 1] . (14) 

This definition leads to the following characterization of 
the boundary of the convex set of A^-outcomes POVM's 

Theorem 4 A POVMP e belongs to the boundary 
of S^N iff o,t least one element Pf of P has a non trivial 
kernel. 

Let's first prove necessity. Consider two different 
POVM's P and Q, and suppose that Ve G [0, 1] P eD 
is still a POVM while P — eD is not. This happens only 
if Ve G [0, 1] Pf — eDf ^ for some /. Then some vector 
■0 must exist such that 

{iP\Pf liP) <e{^\Df\ij), Vee[0,l], (15) 

namely (iplPflip) — 0. Since by hypothesis Pf is positive 
semidefinite, then necessarily ijj e Ker(P/). To prove 
that the condition is also sufficient, consider a POVM el- 
ement Pf with nontrivial kernel, and take e Ker(Py). 
Then consider an event g such that {ip\Pg\il;) > (such 
event must exist for normalization of the POVM), and 
take Df — K\ip){ip\, Dg = —K\i(;){ip\, and — other- 
wise, with K smaller than the minimum eigenvalue of Pg. 
Clearly Ve G [0, 1] P + eT) is a POVM while P-eD is not, 
since the element Pf — eDf is not positive semidefinite. 

We now proceed to study the structure of the faces of 
^N- For such purpose it is convenient to regard a con- 
vex set as a subset of an affine space, whose dimension 



is the number of linearly independent directions along 
which any internal point can be symmetrically shifted. 
Clearly, also the faces of the convex set are themselves 
convex. For example, moving from a point inside of a 
cube one can explore three dimensions while remaining 
inside, whereas, for the cube faces the number of inde- 
pendent symmetric perturbations is two, and for the sides 
this number reduces to one. We will keep in mind the 
above geometrical picture for the classification of the bor- 
der of the convex set of POVM's using the perturbation 
method. 

According to the results of Section |l[ a perturbation 
for a POVM P is a set of Hermitian operators D — {D^} 
with X)e6E^e = 0, and with Supp(Z)e) C Supp(Pe)- Ex- 
pressed in the orthonormal basis of eigenvectors of the 
POVM elements as in Eq. (||), the operators D^. read 

rank(Pe) 
mn—l 

We remind that, according to Theorem ^, non triv- 
ial perturbations for P exist only if the outer prod- 
ucts \vm'^){vn^\ are linearly dependent with e G E and 
1 < n, m < rank(Pe). The total number of such outer 
products is 

r(P) = ^rank(Pe)^. (17) 

eeE 

The number of linearly independent elements in the set 
of outer products \vm^){vn^\ is given by 

/(P) = dim[Span({|z;(„^))(^;(^)|})]. (18) 

Now, as in Theorem |[ we see that the number 6(P) of 
independent perturbations for P is given by 

6(P) = r(P) -;(P). (19) 

In fact, in Eq. (||) we have r(P) = J2eeE rank(Pe)^ vari- 

(e) 

ables -Dnm, whereas the number of linearly independent 
equations is ^(P) = dim(Span({|um^)(i;i'^^|})), whence the 
number of variables which can be written as linear com- 
bination of a linearly independent set is r'(P) — ^(P). 
On the other hand, the dimension of the affinc space 
of the convex set is given by (P{N — 1), since the 
POVM normalization constraint corresponds to (P inde- 
pendent linear equations, with N(P variables. We have 
then proved the following characterization of the border 

of ^AT 

Theorem 5 A POVMV G ii^N belongs to the boundary 
d^N of 3^N iffb{P) < d^{N - 1), 6(P) defined in Eqs. 
(|7^-|7^ being the dimension of the face in which P lies. 

From the above theorem it also follows that a POVM 
P G on the boundary (9^jv of also belongs to 
d3^M with M > N, whence it belongs to the boundary 
of ^. This also implies that d^N C d^M Q d^. 
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V. EXTREMAL INFORMATIONALLY 
COMPLETE POVM'S 

In this section we will give an explicit construction of 
a rank-one informationally complete POVM as in Corol- 
lary ^, in this way also proving the existence of extremal 
POVM's with elements. 

Consider the shift- and-multiply finite group of unitary 
operators 

Upq^Zm\ p,qeZd (20) 

where — {0,1, . . . ,d — 1}, and Z and W are defined 
as follows 

Z^Y.\^®l){j\, W^Y.^^j){j\, (21) 

3 

with denoting the sum modulo d, {|i)}o^^ an orthonor- 
mal basis, uj = , and the sums arc extended to Z^. 
We now prove that the following POVM with out- 
comes is extremal 

Pp,^\Up,vUl, (22) 

for any pure state v = on H satisfying the con- 

straints 

Tv[Ul^v]^0, yp,qeZa. (23) 

In order to prove the statement, first we notice that the 
operators d~^Upq form a complete orthonormal set of 
unitary operators, i. e. they satisfy 

TlppqUp' ^q>] = dSppfSqqi, (24) 

Y^UpqSUlq = dTr[S], (25) 

pq 

for any operator S. From Eq. ( p5| ) it immediately fol- 
lows that Y.pq Ppq = whence {Ppq] is a POVM. Then, 
in order to prove extremality, according to Corollary |^ 
it is sufficient to prove that the d^ operators Ppq are lin- 
early independent, which in turn can be proved by show- 
ing that {Ppq] is itself a complete set in the space of 
operators (completeness along with the fact that the el- 
ements Ppq are implies indeed that they are linearly 
independent). As mentioned after Corollary ^, such a 
kind of completeness for the POVM corresponds to a so- 
called informationally complete measurement^^. The 
completeness of the set {Ppq} is equivalent to the invert- 
ibility of the following operator on H®^ 

F = Y.\Ppq)){{Ppq\, (26) 
pq 

where the double- ket notation is used to remind the 
equivalence between (Hilbert-Schmidt) operators A on H 
and vectors \A)) =A(^ I\I)) of H®^ |/)) e H®^ denoting 
the reference vector |/)) = X^jez^ defined in terms 



of the chosen orthonormal basis By expanding ly = 

IV') ('01 over the basis {Upq} and using the multiplication 
rules of the group, one obtains 

= ^ E e"^^"^ Tr[C/,tHC/.., (27) 

rs 

which allows us to rewrite Eq. ( [26| ) as follows 

F=^Y.\^'^Urs'^^\'\Urs)){{Urs\. (28) 

rs 

Since {Upq))} is an orthogonal basis in H®^, the invertibil- 
ity of F is equivalent to condition (p3|), which is clearly 
satisfied by most density operators = \il))(ip\ (condition 
( p3| ) is satisfied by a set of states \->p) that is dense in H). 
As an example of state satisfying the condition (j2^) , one 
can consider \ip) oc ct^lj) for < a < 1. 

VI. CONCLUSIONS 

In this paper we have completely characterized the con- 
vex set of POVM's with discrete spectrum. Using the 
method of perturbations, we have determined the ex- 
tremal points — the "indecomposable apparatuses" — by 
three alternative characterizations corresponding to The- 
orems and with easier necessary or sufficient condi- 
tions in Corollaries In particular, we have shown 
that for finite dimension d an extremal POVM's can have 
at most d^ outcomes, and an extremal POVM with 
outcomes always exists and is necessarily informationally 
complete. An explicit realization of such extremal infor- 
mationally complete POVM has been given in Section 

The characterization of the convex set of POVM's 
with N outcomes has been obtained by determining its 
boundary d^N, which, in turn, has been characterized in 
terms of the number b(P) of independent perturbations 
for the POVM P in Eq. (|l^). This has lead to a sim- 
ple characterization of the boundary in terms of simple 
algebraic properties of a POVM lying on it. Since 
is also a subset of the boundary of the full convex 
set ^ of POVM's with discrete spectrum, our result also 
provides a complete characterization of 

Finally, in Appendix ^ we reported an algorithm for 
decomposing a point in a convex set into a minimum 
number of extremal elements, specializing the algorithm 
to the case of the convex set of POVM's with N 
outcomes. 



APPENDIX A: ALGORITHM FOR 
DECOMPOSITION OF INTERNAL POINTS 
INTO EXTREMALS 

In this appendix we provide an algorithm to decompose 
a POVM P e into extremal ones. We first present 
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the algorithm in the general case of an abstract convex 
set, and then we specialize it to the case of POVM's. 

Consider a convex set and a point p inside it. We 
want to decompose p into extremal points. We first need 
two ingredients, which depend on the specific convex set 

under consideration: the affine space si in which 
is embedded (this is just the space of legitimate "per- 
turbations"), and an "indicator" which is positive 
for p e zero on the boundary, and negative outside 
. Now, starting from p inside '£ ^ we move in some 
direction d in jz/ until a face of is encountered at 
A-(- = max{A : p + Ad G "^j (using our indicator A4- 
is given by the value of A where i(p + Ad) changes sign). 
Similarly in the opposite direction one hits the boundary 
at A_ ^ max{A : p — Xd G "^}. The point p can now be 
split into the convex combination 



A_ 



A, 



P 



X+ + X- 



A+ + A_ 



p± = p ± X±d. (Al) 



If p was in the interior of any face of the boundary can 
be encountered, while if p was already on the boundary 
of the perturbation d brings p on a " face of a face" , 
e. g. moving on a face of a cube toward an edge. In any 
case, the dimension b{p) of the face to which p belongs is 
decreased at least by one. 

By applying the same splitting scheme to both p± re- 
cursively, we obtain a weighted binary "tree" of points 
rooted in p, with the property that the point p' at each 
node can be written as convex combination of its descen- 
dants, and with a depth bounded by the dimension b(j>) 
of the face to which p belongs. Of course the "leaves" 
of the tree are extremal points pi , and one can combine 
them to obtain the original point as p = (XiPi weight- 
ing each leaf pi with the product of all weights found 
along the path from the root p to the leaf p,;. Unfor- 
tunately, this raw algorithm can produce up to 2^ ex- 
tremal points Pi, for dimension r of the affine space of 
each leaf being addressed by the vector di — p — pi, with 
y^aidi = 0. However, by the Caratheodory's theorem 
||16| , we know that at most r+1 extremal point are needed 
to decompose p. Indeed, if the number of di is larger than 
r + 1, then they must be linearly dependent, and there 
must exist Xi 's not all vanishing and not all positive such 
that = 0. Since J^t^^A ~ PiJ^i^A) = 0; by 

choosing the greatest ^ such that a.i — ^Xi > Vi, one 
finds that can be written as a convex combination of 
a smaller number of d^'s. This procedure can be applied 
repeatedly to the remaining d^'s (the aiS must also be 
upgraded) until their only combination giving is the 
one whose coefficients are all positive: at this point their 
number is for sure not larger than r + 1. Therefore, from 
an intial decomposition with many elements, we end up 
with a decomposition of p into at most r -\- 1 vectors di 
and probabilities Ui such that '^i^i ctidi = 0, whence 

p = J2l=i ^iip + di) = OLiPi. Notice that the eval- 

uation of A± at each step involves an eigenvalue eval- 
uation, whence the algorithm generally doesn't provide 
analytical decompositions. 



In order to specialize the algorithm to the case of the 
POVM's convex set 3^Ni we need to specify both the 
corresponding affine space jz/at and the indicator l of the 
border. The affine space is the real d^(7V— l)-dimensional 
linear space of vectors D — {D,.} of N Hermitian opera- 
tors De, with J2eeE = 0. This can be obtained as the 
real span of projectors \n, e){n, e\ along with ^\n, e)(m, e| 
and e)(m, e| for n = l,...,d and n < m, where 
X — '^RX + i'<sX is the Cartesian decomposition of the 
operator X, {|rt, e)} denotes any orthonormal basis for 
the e-th copy of the Hilbert space H of the quantum 
system with d = dim(H), and e G E' = E/{1}, and 
Di — — ^ggg, diagZ^e- However, if such global basis 
is used, then when the search algorithm starts from a 
POVM which is already on a face of the convex one has 
the problem that generally the basis is not aligned with 
the face itself, and for a generic direction the perturbed 
POVM either exits from the convex, or it moves inside 
it. For that reason it is convenient to consider a "lo- 
cal basis" of perturbations for a given P. This can be 
constructed by considering the set {Xm\i} of Hermitian 
operators defined as follows 



X 



(e) 



•J V, 



n < m 



(A2) 



Then pick up ?(P) linearly independent elements, which 
we will denote by Vmr. 
so that we can write 



we will denote by Vmn , and call the remaining ones Wmn , 



^ ^ ran / j ^epQ pq — / ^ / ^ epq pq ' \ ) 

e^f PI 



epq 



where the second identity is a consequence of linear in- 
dependence of operators {\vrn){vn^\} for fixed e. Then 
we can construct the following basis {De{fmn)} for P 
perturbations 



Djifmn) = Y,cj;;v(p-Wifl 



pq 



(A4) 



Clearly one has X^ees Deifmn) — yfmn, and modulo 
a suitable rescaling one has P^. ± De{fmn) > \l fmn. 
Notice that, by construction, the operator s D ejfnm) are 
linearly independent. In fact, using Eqs. (A4), a generic 
linear combination of De{fmn) for each fixed e S E will 



result in a linear combination of {Fmn} U {W^l''} for the 

same e. But the set {Fmn} U {Wjm^ is linearly indepen- 
dent for fixed e G E, due to Eq. (|A4|). 

As regards the indicator of the boundary, this is sim- 
ply the minimum of the eigenvalues of the operators Pe, 
which changes from positive to negative when P crosses 
the border. 
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